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ABSTRACT 

The dynamical evolution of planetary systems leaves observable signatures in 
debris disks. Optical images trace micron-sized grains, which are strongly affected 
by stellar radiation and need not coincide with their parent body population. Ob- 
servations of mm-size grains accurately trace parent bodies, but previous images 
lack the resolution and sensitivity needed to characterize the ring's morphology. 
Here we present ALMA 350 GHz observations of the Fomalhaut debris ring. These 
observations demonstrate that the parent body population is 13-19 AU wide with 
a sharp inner and outer boundary. We discuss three possible origins for the ring, 
and suggest that debris confined by shepherd planets is the most consistent with the 
ring's morphology. 

Subject headings: Planet-disk interactions — planetary systems — Submillimeter: 
planetary systems 



1. Introduction 



The Fomalhaut debris system is a natural laboratory for testing planet formation theories. 
The nearby (7.69 pc. Ferryman et al. 1997) A3V star (Di Folco et al. 2004) is surrounded by 
an eccentric debris ring with a peak brightness in scattered optical light at a semi-major axis 
a ~ 140 AU ( Kalas et"aL]|2005 i henceforth K05). The inner edge is sharply truncated, which, 
along with the ring's eccentricity, suggests that a planet is shaping the ring's morphology ( |Wyatt| 
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eraL][T999t |Quillen|[2006t [Chiang et aL][2009] ). A candidate for the reputed planet (Fom b) has 
been discovered in the optical ( [Kalas et al.|2008[ ). The observations are consistent with a super- 
Earth mass planet embedded in a planetesimal swarm ( Kennedy & Wyatt|201 1 ). 



Independent of direct detection of any planet, observations of the ring's morphology can 
constrain the dynamical history of the Fomalhaut system, including properties of any planets 
near the ring. The effect of radiation pressure on a dust grain's orbital eccentricity, assuming 
an initially circular orbit, is e = (2pss/(pooSoo) — 1) \ where ps is the grain's internal density 
and s the grain radius. Here, oo represents grains that are unbound by radiation pressure. For 
Fomalhaut, Soo ~ 8yum for poo ~ 1 g/cc. Most of the scattered optical light should be from 
the smallest bound grains (between ~ 8 and 16pm, see [Chiang et al.|2009[ ). The expected free 
eccentricity of a 16pm grain due to radiation pressure is ~ 0.3 (for p^ = poo), while grain sizes 
near 1 mm will have free eccentricities e < 0.01, making mm grains excellent tracers of parent 
bodies. While previous observations of mm grain emission do detect large grains, they lack 



the resolution needed to characterize the parent body morphology (Holland et al. 1998, Ricci 
et al.|2012 ). Here, we present high-resolution 850 pm ALMA images that resolve the northern 
section of Fomalhaut's ring. 



2. Observations and Reduction: 

Fomalhaut's ring was observed using ALMA cycle in the compact configuration, measur- 
ing projected baselines from 14 to 175m (Table[T]). Observations were centered on the expected 
position of Fom b at RA = 22h:57m:38.65s and 5 = -29d:37':12.6" (J2000, proper motion in- 
cluded). The total on-source integration time was 140 min. The observations were taken at 357 
and 345 GHz (in the upper and lower sidebands) using the Frequency Domain Mode in dual 
polarization with 4x1875 MHz bandpasses (2 in each of the sidebands). Neptune was used 
as an absolute flux calibrator, and J 1924-292 for bandpass calibration. Atmospheric variations 
at each antenna were monitored continuously using Water Vapor Radiometers (WVR), as well 
as regular hot/ambient load measurements. For time-dependent gain calibration, the nearby 
quasar J2258-279 was observed every 8 minutes. Data were reduced using CASA 3.4 ( Mc-| 



MuUin et al. 2007): calibration involved removing the effects of rapid atmospheric variations 
at each antenna using the WVR data, correcting the time- and frequency-dependence of system 
temperature, and correcting for the complex antenna-based bandpass and time-dependent gain. 
Amplitude calibration used the CASA Butler- JPL-Horizons 2010 model for Neptune, which 
gives an estimated systematic flux uncertainty of ±10%. The calibrated measurement set was 
spectrally binned to a channel spacing of 49MHz, and then CLEANed using the Cotton-Schwab 
algorithm, combining all channels to give the final continuum image. The primary beam cor- 
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rection (PBC) was performed using the voltage pattern for an Airy function in with an effective 
dish diameter of 10.7 m, a blockage diameter of 0.75 m, a maximum radius of 1.784°, and a 
reference frequency of 1 GHz (see CASA functions vp.setpbairy and sm.setvp). 



3. Results 

In Figure 1 we present the cleaned and primary beam-corrected images. We refer to these 
as the uncorrected and corrected images, respectively. In the uncorrected image, the peak bright- 
ness is ~ 0.84 mJy beam~^, the total flux density of the image is 20.5 mJy, and the RMS is 
~ 60 /iJy beam~^ (0.32 mK|^ The synthesized beam is ~ 1.5" x 1.2" with position angle 77°. 

The corrected image shows all emission to 7% power, which includes the ring and all 
apparent emission related to the star. The total flux in the corrected image is 45.5 mJy. The 
projected image accounts for approximately half of the ring, so assuming equal brightness for 
the second half and accounting for the expected contribution from the star (3 mJy), we find a 
total flux density of ~ 85 mJy, reasonably consistent with previous measurements (81 ± 7.2 and 
97 ± 5 mJy [Holland et al.|1998[|2003| ). Almost all of the emission is confined to a thin ring. 

The ring's surface brightness reaches a maximum of 1.54 mJy beam~^ near the ani'fil^and 
remains bright through the inferred apocenter. There also appears to be a reduction in the ring 
brightness to the SW. We will discuss the significance of these azimuthal variations below. The 
star's peak brightness is 3.4 mJy beam"^ and has a flux density ~ 4.4 mJy (3 mJy expected 
for a blackbody photosphere), which suggests possible excess emission. We caution that the 
measurement of total stellar flux is strongly dependent on the PBC. Nonetheless, a detection of 
excess emission is consistent with the results of ( Stapelfeldt et al.|2004 ). The PBC error would 



need to be 47% at the 10% power point to account for the excess, more than 4 times the ALMA 
beam specification which has yet to be completely verified. 

Without the second ansa, we cannot reliably determine the ellipse. We show in Figure 1 
two ellipse segments that are based on fits to optical images (K05).The center for both ellipses, 
denoted by the plus sign, is 0.29" W and 1.7" N of the star. The ellipses have a semi-major axis 
a = 18.31" and a position angle (PA = 336°). Assuming the ellipses represent sky-projected 
circles, the red ellipse has an inclination of i = 66° (K05) and the blue 67°, which more closely 



'Unless otherwise stated, relative errors are assumed to be equal to the RMS value of 60^ Jy beam ^. Flux 
densities are only expected to be good to about 10%. 

^Ansa refers to the sky-projected section of an astrophysical ring that is farthest from its central body, tradia- 
tionally referring to planetary rings. 
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corresponds to our best-fit sky models (below). 

Surface brightness profiles of the ring (Fig. 2) are taken along the 7 shces shown in Figure 
1 using bilinear interpolation. Each slice can be fit by a single Gaussian except slice 4, which 
requires a double Gaussian due to elongated emission. Altogether the slices, excluding slice 4, 
have a combined FWHM ~ 2.19", corresponding to a projected parent body width ~ 17 AU, 
well-resolved by the synthesized beam. 

Deprojected slices through the corrected and uncorrected images at the ring's ansa are 
shown in Figure 2, which includes all points within ±2° of the PA. The FWHM = 1.87 ± 0.03" 
and 1 .85 ±0.03" for the corrected and uncorrected profiles, respectively. The surface brightness 
in the corrected image peaks at ~ 18.4" (141.5 AU). The peak is similar to that found for the 
small grains seen in the optical light radial brightness profile (K05), but the mm grains are much 
more tightly confined, consistent with radiation pressure effects. 

We estimate the mm grain mass by assuming the ring is optically thin and comprised of 
grains with s ~ 1 mm and ps ~ 2.5 g/cc. The grain temperature Tg ^ 0.7Tstar(-Rstar/-D)^/^, 
where Tgtar = 8750 K and -Rstar = 1.82i?0 are the stellar temperature and radius ( Di Folco] 



eraLl[2004] ), respectively. At D = 140 AU, Tg ^ 48 K. The ring's total flux density of 82 mJy 



(excluding the expected 3 mJy from the star), so ~ 0.017 M©. 

If collisions are the main removal mechanism of mm grains, then a lower limit can be 
placed on the parent body mass Mpbdy/^age ~ ^mm/^cou for ring age tagc and mm-grain colli- 
sion timescale tcou ~ (rfi)~^. The vertical optical depth through the ring can be estimated by 
r ~ (surface brightness/flux density per grain) (vrs^). The peak of the ring's surface bright- 
ness is ~ 1 mJy beam^^. The synthesized beam effective area is ~ 2 AU^ at 7.69 pc. For 
s ~ 1 mm and Tg ~ ASK, r ~ 10"^. At 140 AU, tcou ~ 2 Myr. Combined with the above 
mass estimate, the ring should be losing mass at a rate of M ~ O.OlMEarth Myr^^, requiring 
Mpbdy > 1.7MEarth (for tage = 200 Myr), consistent with ( [Chiang et al.|2009[ ). 



Best-fit models for the spatial distribution of mm grains are selected by minimizing the 
difference between model and ALMA visibilities, where the u-v plane sampling from the ob- 
servations is used in CAS A to observe each model (tool sm.predict). Ring models are produced 
using the following assumptions: (1) The ring is circular and has an offset relative to the star 
that is in agreement with the deprojected ellipse models (see Fig. 1), approximating a grain 
orbital eccentricity ~ 0.1. (2) The ring is uniform in azimuth. (3) The radial profile for mm 
grains is described by either a Gaussian distribution in semi-major axes or a power law in sur- 
face density. For the Gaussian, the FWHM is varied from 9.4 to 21 AU, and the inner radial 
half-maximum point from 131 to 140 AU. For power law surface density profiles, the annular 
width is varied from 10 to 30 AU, the inner edge from 133 to 138 AU, and the power law index 
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from -9.5 to -1.5. (4) The vertical distribution for grains follows an exponential decay with a 
scale height given by a constant angle above the midplane between and 2.5°. In each model, 
the emission from grains at a given distance from the star is calculated from the Planck function 
with the grain temperature Tg determined as described above. The total flux density for the ring 
at 850;um is normalized to either 80, 90, or 100 mJy. The star's flux is set to 3 mJy. Each model 
is projected onto the sky with PA = 336° and i = 66 to 67.5°. The model ring's apocenter is 
set to the inferred apocenter of Fomalhaut's ring. It should be noted that models of a uniform 
ring can become brighter near the ansae due to projection effects. Smooth, symmetric models 
cannot, however, account for the excess brightness near apocenter. Such brightening could be 
explained by variations in the mass distribution of grains or by variations in ring's radial thick- 
ness. Collisions between parent bodies could also produce asymmetries ( [Wyatt^ Dent||2002| ). 



The best-fit models with a Gaussian semi-major axis distribution have a scale height cor- 
responding to an opening angle of 1.0° ± 0.25, a FWHM = 16 ± 3 AU, and an inner radial 
half-maximum at 13515^ ^U. The 1-a errors for model parameters are estimated by including 
models for which Ax^ < 1. Models with i = 66.75° and a total flux density of 80 mJy are 
preferred. Most of the dimming in the SW is consistent with a loss of sensitivity (Fig. 1). The 
appearance of a hole may be due to brightening in the southern-most part of the ring, which 
may be due to noise. The NE is not subtracted in the residual image (Fig. 1, bottom-right) and 
contains contiguous excess surface brightness. The star is also not well- subtracted, suggesting 
possible excess. Some of the extended surface brightness noted in slice 4 of Fig. 1 remains, but 
is consistent with a high-noise peak. 

Some select, steep power law surface density profiles are also consistent with the data. 
These are limited to a surface density oc r~^-^, an inner edge of 135^| 5 AU, an opening angle 

I. 0° ± 0.5, and a preferred i = 67.25°. For these power laws, the half-maximum width is about 

II. 4 AU. 



4. Discussion 

Any formation mechanism for Fomalhaut's ring must address the following constraints: 
(a) The ring has an eccentricity ~ 0.1. (b) The vertical scale height of the parent bodies is 
consistent with an opening angle of ~ 1.0° from the midplane. (c) The FWHM of the parent 
bodies for the Gaussian model is ~ 16 AU, giving a width-to-height aspect ratio ~ 7. (d) The 
outer edge of the parent bodies is consistent with being as sharply truncated as the inner edge, 
(e) The micron grains are more radially extended than the parent bodies, but this is satisfied for 
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a variety of models due to radiation pressure. We consider three formation models. One is that 
the ring is sculpted by a single interior planet ( |Wyatt et al.|1999^[Kalas et al.|2005^[Quillen|2006 



[Kalas et al.|2008^ [Chiang et aL]|2009| ). To account for requirement (d), this solution requires an 



abrupt truncation of the outer disk through, for example, a stellar flyby ( [Ida et al^|2000p . We 
will return to this scenario below. 

Another possible model is that the ring is a remnant from a collision between two plan- 



ets (Mamajek & Meyer 2007). Unless constrained, the radial width of the ring will slowly 
spread due to collisions between ring particles with different semi-major axes. Take the ver- 
tical velocity dispersion 5vz ~ vxsmi for Keplerian orbital speed vk- If the radial velocity 
dispersion is of the same order, then the typical a difference between two colliding particles is 
6a ~ aSe to factors of order unity, which for free eccentricity Se ~ sin i is the scale height of 
the disk, h. The effective coUisional viscosity of the ring Uc ~ /i^/tcoii- Using a scale height 
of 2.5 AU and a coUisional time between mm grains of 2 Myr, Uc ~ 2.2 x 10^'^ cm^/s, which 
means the ring is diffusing on timescales ~ 90 Myr. If the ring is younger than the star, no 
planetary influence is necessary apart from the dynamics that led to the collision. However, 
shattering gravitationally bound objects requires impact speeds several times the mutual escape 
speed = 2G(Mi + M2) / {Ri + R2), for bodies with masses Mi and M2 and radii Ri and 
R2). ( |Asphaug|20T0t|Leinhardt & Stewar"t||2012| ). For two IMEarth planets (Afpbdy > IMEarth 
for tage ~ 90 Myr), Vesc ~ 11 km/s. The orbital escape speed at 140 AU from Fomalhaut is only 
about 5 km/s. The impact speeds necessary to destroy Earths are not attainable. 

Finally, we propose a third formation model in which the ring's morphology is dominated 
by shepherd planets, analogous to the shepherd moons Cordelia and Ophelia of Uranus's e 



ring (Goldreich & Tremaine 1979, Smith et al. 1986). Here, the ring is confined by angular 
momentum exchange between ring particles and the planets. We explore this possibility with 
N-body simulations, using the Bulirsch-Stoer algorithm in Mercury ( Chambers]| 1999 ). In each 
simulation, two planets of mass Manner and Mouter are placed at amner and aouter relative to 
the 2.3 Mq star. Their relative inclinations are set to zero. Planet eccentricities are set to force 
the planetesimals (see method in Chiang et al. 2009) at ~ 143 AU to have the eccentricity 
observed for the ring (~ 0.11). For each run, 10^ massless particles are placed between Tinner 
and router with a surface density profile oc r"^, where Tinner = max(125 AU, ainner + '^ AU) and 
router = min(160 AU, aouter — 1 AU). The integrations run for 10* years. Massless particles 
are removed from simulations if the particle collides with a planet or if its radius from the star 
exceeds 10^ AU. 

The masses of the planets set the steepness and skew of the ring's radial profile (Fig. 3). 
Mutual perturbations will cause the system to evolve, which narrows parameter space by de- 
manding that interactions between the planets do not destroy the ring. In the massless ring- 
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particle limit, either both planets must be < SMgarthj or they must have an extreme mass ratio. 
The planets do not need to be in resonance. A super-Earth and a Mars-mass planet produce 
the most narrowly-peaked Gaussian for the radial parent body distribution among the simula- 
tions presented here. As the masses for shepherd candidates are comparable to the minimum 
estimated parent body mass in the ring, self-gravity of the ring may play a role in the system's 
evolution. If the scattered optical light observations have detected Fom b, then the candidate 
could be the innermost shepherd, although there is a discrepancy of ~ 10 AU between Fom b's 
proposed a and the inner shepherd in our models. 

Fomalhaut's debris has some similarities to the Kuiper belt. Both systems are thin rings, 
Ar/a ~ 0.1, which seem to owe their morphologies, at least in part, to the presence of planets. 
There are nonetheless substantial differences. The outer edge of the Kuiper belt is near the 2:1 
resonance with Neptune (Trujillo & Brown 200 For the outer edge of Fomalhaut's ring to 
be at the 2:1 resonance with an interior planet, the planet would need to have an a ~ 95 AU. 
To truncate the ring at this location, the planet would need to be > 3Mj ( Chiang et al.||2009 ) 
and could have been decected by Spitzer ( [Janson et al.||2012] ). If the ring had formed due to 
outward migration of an interior planet, we would expect particles to be trapped in resonances, 
like the Plutinos in the 3:2 resonance with Neptune. We detect only a single ring. The Kuiper 
belt has a vertical thickness given by a dynamically hot and cold population with Gaussian 
widths a 



17° ±3 and a 



2.2° +0 6, respectively (|Brown|200l|). While the Kuiper beh object 



inclination distribution seems to be consistent with truncation due to a stellar flyby (Ida et al. 



2000), Fomalhaut's parent body population may be too cold to be explained by this mechanism. 
This leads us to favor shepherd planets as the explanation for the ring's morphology. 
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Table 1: Observing Log: Target Fomalhaut's Ring and Fomalhaut b. The precipitable water 
vapor (PWV) is given in mm. 



Date 

dd/mm/yyyy 


Time 
(UTC) 


On-Source 
(min) 


Antennas 


PWV 


22/09/2011 


23:45:00.0 


25 


13 


0.63 


23/09/2011 


00:49:30.0 


25 


13 


0.48 


23/09/2011 


02:08:15.3 


25 


13 


0.80 


18/10/2011 


03:31:00.9 


25 


15 


0.62 


18/10/2011 


04:42:55.6 


20 


15 


0.63 


18/10/2011 


05:44:40.4 


20 


15 


0.57 
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Offset West (arcsec) 




Offset West (arcsec) 



Fig. 1. — Top-left: 350 GHz ALMA image of Fomalhaut's ring. The RMS noise is 
~ 60 /iJy beam~^, and the contours represent 0.24, 0.33, 0.47, 0.66, 0.92, and 1.4mJy beam^^. 
The circles are centered on the pointing center, and their diameters show the beam's half and 
7% power. Ellipses with PA = 336°, and i = 66° (red) and 67° (blue) are shown, with the 
ellipse centers given by the plus sign. The slices labeled 1-7 are used to show surface brightness 
profiles in Figure 2. The coordinate axis is centered on the star, which has a peak brightness of 
0.49 mJy beam"^ (uncorrected). Top-right: The primary beam corrected image. The ring be- 
comes bright near the ansa, and remains bright through the inferred apocenter. The contours are 
the same as in the top-left, but begin at 0.47 mJy beam"^ Bottom-left: The best-fit skymodel 
is compared with the dirty ALMA image. The grayscale and contours show mJy beam"^, with 
the contour levels labeled on the colorbar. Bottom-right: The residual of the data minus the 
model. The model subtracts the emission well, but there remains excess emission in the NE and 
around the star. 
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AU From Gaussian Peak AU From Ellipse Center 
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Arsec From Gaussian Peak Arcsec From Ellipse Center 

Fig. 2. — Left: Normalized surface brightness profiles for slices 1-7 in Fig. 1. All profiles 
except slice 4 can be fit by a single Gaussian. A fit through all slices, excluding slice 4, is 
shown with the dot-dot curve. Right: The deprojected radial surface brightness profile through 
the ansa for the corrected and uncorrected data. Black curves show corresponding Gaussian 
fits. The heavy-dashed curve shows the scaled scattered optical light profile (K05). The dot- 
dot curve shows the scaled Gaussian fit for the combined slices in the left panel. 
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Fig. 3. — Radial profiles for a ring of parent bodies between two planets. Left-to-right: 



Mi, 



QM Earth and Mi, 



^M Earth- Top-to-bottom: Mouter = Manner, M, 



outer 



Minner/^, and Mouter = firmer /lO. Each color represents a profile-planet pair. The gray 
dot-dashed line represents the best fit to the observed radial profile. The left y-axis shows the 
eccentricities of the planets, and the right shows the number surface density of particles, with 
each curve's peak normalized to unity. The range of eccentricities occupied by the planets is 
indicated by the vertical error bars. The x-axis shows the semi-major axes for the planet and the 
radial distance of the planetesimals from the ellipse center. 



